SOME APPLICATIONS OF SEMI-DISCRETE VARIATIONAL 
INTEGRATORS TO CLASSICAL FIELD THEORIES 
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Abstract. We develop a semi-discrete version of discrete variational mechanics with applica- 
tions to numerical integration of classical field theories. The geometric preservation properties 
are studied. 



1. Introduction 

The calculus of variations is a fundamental tool in the description and understanding of Clas- 
sical Mechanics and it is an area of active research. One part of this activity was dedicated 
to uncover the geometrical structures behind such formalisms. Many physical systems not only 
evolve in time, as in Classical Mechanics, but also posses a continuous spatial structure. This is 
the setting of Classical Field Theories in both the Lagrangian and Hamiltonian formalisms [H [3] . 

One traditional way to analyze these problems has been to pass to the Hamiltonian formalism 
(or Lagrangian formalism) using a space-time decomposition of the parameter space, and then 
applying classical Dirac's theory of constraints (in the singular case, which is typical in field 
theories). In this way, we obtain a well-posed system of equations of motion that can be even- 
tually integrated or numerically simulated. It is clear that using this space-time decomposition 
we broke the original covariance of the theory and perhaps some geometrical structure is lost, 
but the treatment of the equations is, in many aspects, more easy and some of the geometri- 
cal structure is still preserved. Moreover, for numerical simulation of the equations of motion, 
after the space-time decomposition, we eventually obtain a Hamiltonian system and symplectic 
integration methods may be useful to solve numerically the initial problem. 

In this sense, it may be useful to introduce geometric integrators, that is, numerical schemes 
which preserve some of the extra features of geometric nature of the dynamical systems. Usually, 
these integrators can run, in simulations, for long time with lower spurious effects (for instance, 
bad energy behavior for conservative systems) than the traditional ones [SI E] ■ 

A particular case of geometric integrators are variational integrators. These integrators have 
their roots in the optimal control literature in the 1960's and they enter in the "geometric 
differential arena" after the pioneering work of Veselov [19] and Moser and Veselov [TSl US] • In 
these papers, there appears the discrete action sum, discrete Euler-Lagrange equations, discrete 
Noether theorem... These integrators have been adapted for Field Theories pj [10]. All these 
integrators have demonstrated exceptionally good longtime behavior. 

In this paper, we will develop the theory of semi-discrete variational integrators for Classical 
Field Theories. The basic idea is to consider an spatial truncation that reduces the partial differ- 
ential equations derived from the Euler-Lagrange equations to a system of ordinary differential 
equations [HQS]. The main objective of the paper is to study the geometric properties after this 
spatial truncation: preservation of forms, energy preservation and momentum preservation. 
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2. Discrete variational calculus 

First, we will recall discrete variational calculus, following the approach in |12| and references 
therein. A discrete Lagrangian is a map L d : Q X Q — > M, which may be considered as an 
approximation of a continuous Lagrangian L : TQ — > M. Define the action sum S d : Q N+1 — > K 
corresponding to the Lagrangian by 

iv 

S d (q ,...,q N ) =22L d (q k -i,qk) , 

k=l 

where % G Q for < fc < N. 

Observe that for any covector a G X2 \{Q x Q), we have the decomposition a = oil + oli 
where on G T*.Q, thus, 

dL d (q , q{) = Di_L d (q Q , g x ) + D 2 L d (q ,qi) . 

The discrete variational principle states that the solutions of the discrete system determined 
by must extremize the action sum given fixed points qo and gj\r. Extremizing over gfc, 
1 < k < N — 1, we obtain the following system of difference equations 

(2.1) DiL d (q kl q k+1 ) + D 2 L d (q k -i,q k ) = . 

These equations are usually called the discrete Euler- Lagrange equations. Under some regularity 
hypothesis (the matrix (Di2L d (q k , q k +i)) is regular), it is possible to define a (local) discrete flow 
T : Q x Q — ► Q x Q, by T(gfc_i, %) = (q k ,qk+i) from Equations (|2"7T]) . 
Define the discrete Legendre transformations associated to L d by 

FL- : Q x Q — > T*Q 

(«o,«i) 1 — > (go ! --Diid(g , o,gi)) , 

FL+ : Q x Q — > T*Q 

(?o,9i) 1 — »■ (ffij^idCgoffi)) > 

and the 2-form = (FL^)*ujq = (FL^)*ojq, where is the canonical symplectic form on T*Q. 
The discrete algorithm determined by T preserves the symplectic form uj d , i.e., T*^ = u> d . More- 
over, if the discrete Lagrangian is invariant under the diagonal action of a Lie group G, then the 
discrete momentum map J d : QxQ -> g* defined by (J d (qk,%+i),0 = {D2L d {q k ,q k+ i),^ Q (q k+1 )) 
is preserved by the discrete flow. Therefore, these integrators are symplectic-momentum preserv- 
ing integrators. Here, £q denotes the fundamental vector field determined by £ G Q, where q the 
Lie algebra of G. 



3. Variational Calculus in Classical Field Theories 

Consider a locally trivial fibration tt : Y — > X, where Y is an (to + n + l)-dimensional 
manifold and X is the parameter space usually equipped with a global decomposition in space 
and time, that is, X — R x P, dimP = n. We shall also fix a volume form on X, that will be 
denoted by to We can choose fibred coordinates {x^,y l ) in Y, so that ■K(x fJ ',y l ) = (x^), where 
(x M ) = (x°, x 1 , . . . , x n ) and x a = t represents the time evolution. Assume that the volume form 
is r] = dx° A ... A dx n . Here, < fi, i/, ... < n and 1 < i, j, ... < to. We shall also use the following 
useful notation d n x^ := io/dx^V- 

The first order jet prolongation J 1 tt is the manifold of classes j^<f> of sections <fi of tt (whose set 
will be denoted by r(7r)) around a point x of X which have the same Taylor expansion up to order 
one. J 1 ^ can be viewed as the generalization of the phase space of the velocities for Classical 
Mechanics. Therefore, J l ir, which we shall denote by Z ', is an (n+l +m+ (n+ l)TO,)-dimensional 
manifold. If we have adapted coordinates (ir M , y l ) in Y, then we have induced coordinates in Z , 
given by 

/ [j. i i \ /il n i i i i \ 

\ x j V j z fj.) — l 1 ! * , y , Z t , Zi, . . . , Z n ). 
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Suppose that we are given a function L : Z — > R of class C 2 in its (n + 1 + m + (n + l)m)- 
arguments. Let be a compact (n+ l)-dimensional submanifold on X. We can thus construct 
the following functional 

3l(4>) = [ L(x»,y l (x),yl(x)) dx° A ... A dx n 

JG X 

for all section cj> G r(7r), where (/>(x) = y l (x), y l u(x)) 

Definition 3.1. A section <f> G is a solution of the variational problem determined by L if 
and only if <j> is a critical point of 3l- 

Extremizing the functional 3l we obtain the Euler-lagrange equations: 

dL d f dL\_ {] 0<i<m 
By 1 dx^ \dzij 

Now, denote by := J^- and by p := L — ztp^, then for a given Lagrangian function L and 

a volume form r\ we can construct the Poincare- Cartan (n + l)-form 

/ BT \ BL 

B L = (L-z l —)d n+1 x + —dy l Ad n x^ 

= (pdx» +ffidy i ) Ad n x" 

From this form, we can also define the Poincare- Cartan (n + 2)-form as £Il '■= — dO^. In 
induced coordinates is expressed as follows 

n L = —(dp A dx^ + dp^ A dy*) A d n x fl 

^From these equations it is easy to derive an intrinsic version of Euler- Lagrange equa- 
tions. In fact, a section <\> G r(7r) is an extremal of 3l if and only if 

(iVra e o i )=o 

for every vector field £ on Z. 

We refer the reader to reference [2] for more details and a complete derivation of the equations 
and the geometric framework for classical field theories. 

4. Motivating Example: The nonlinear wave equation 

Consider the nonlinear wave equation given by 
(4.1) Uh = S x a'(u x ) - /'(«) 

where u : U C R 2 — > M and a, f are smooth functions. If a(u x ) = u 2 /2 then we obtain the 
semi-linear wave equation. 

Equation (|4.1|) corresponds to the Euler-Lagrange equation for the lagrangian function 

L(u,u u u x ) = \^u\ - a(u x ) - f(u) 

Observe that, in this particular case, the Lagrangian does not depend on the parameter space 
I = lxM. 

Now, replace the x-derivative in the Lagrangian by a simple difference (for simplicity, we will 
work with a uniform grid of TV + 1 points, h = L/N) as follows: 

T i r. \ r. \ \ 1 f Mt + (ui)t \ 2 ( ui-ug \ ( u x + u \ 
L ad {uQ,ut, {u )t, \ui)t) = - I ^ ) ~° 1 — h — ) ~J { — 2 — J 

In this case, the Lagrangian is a function L s( i : TR x TR ~ T(R x R) — ► R. This spatial- 
discretization is useful, since as we will show in the next section, applying a suitable variational 
principle we obtain a semi-discretization of Equation (|4.ip , replacing the partial differential equa- 
tion by a system of ordinary differential equations. After this step, the equations of motion can 
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be integrated in time or numerically integrated using, for instance, a symplectic method as the 
symplectic Euler scheme. 

In the next sections, we will show that with these semi-discretizations some of the geometric 
properties are preserved. 

5. Semi-discrete variational calculus 
Given a smooth manifold Q consider the following sets 

£n,[o,t] = {y:{0,...,iV}x[0,T]^g|y(fc ! .)eC 2 ([0,T])forallfc 

and the curves t — > y(0, t) and t — > y(N, t) are fixed} , 
e N ,[o,T] = {y.{0,...,N}x[0,T}^Q\y(k,-)eC 2 ([0,T}) 

and the values y(k,0) and y(k,T) are fixed for all k} . 

The choice of one of these sets depends of the different boundary conditions of the initial problem. 
Of course, another boundary conditions can be analyzed considering suitable adaptations of both 
situations. 

Now, let L : TQ x TQ ~ T(Q x Q) -> M be a Lagrangian function. If (« ,«i) G T 9o Q x T gi Q 
then, as in Section [21 we have that 

dL(v ,vi) = (DiL)(v ,vi) + (D 2 L)(vo,Vi), 

with (DiX)(«o,«i) G T; o (TQ) and (P a i)(«0,«i) G T* (TQ). 

On the other hand, if (g§) (respectively, (g|)) are local coordinates on an open subset Uq 
(respectively, (t/i)) of Q such that go G Uq (respectively, gi G C/i), then we may consider the 
corresponding local coordinates ((/q, gg, q\, q\) on TQ x TQ and it follows that 

T)iL(u ,wi) = ^2(ttj dqliivo^i) + — dq^VQ^^), 

i=i ^ok^i) ^ok^o^!) 

- dL dL 
D 2 L(v Q ,v 1 ) = ^2(ttj dq\{v ,vi) + — dq\(v ,Vi)). 

Moreover, we will use the following notation 



DiL(v ,vi) 


E dL 


rf9o( u o,wi) 


D 2 L(v Q ,v 1 ) 


i=l ^0|(« ,«i 




D 3 L(v ,vi) 


Am 








D 4 L(v ,vi) 


E 5i 


dgi(uo,«i) 



Note that DiL(vo,v\) and D^L{vq,v{) may be considered as 1-forms on Q at the points go 
and gi, respectively. In fact, 

(M*°) 4 (E7TT d «o(9o)) = D 1 L{v ,v 1 ), 
i=l °«0 |(* ,t<i) 

((^)* 1 )*(E7TT <M(«0) = «3i(«o,«i), 

i=l ^1 l(«o,fi) 

where ((r Q )^)* : T g *Q -> T* o (TQ) (respectively, ((tq)?)* : T* Q -» T^(TQ)) is the dual map 
of the linear epimorphism {tq)Z° : T VQ (TQ) -» T go Q (respectively, (tq)^ 1 : T Vl (TQ) -> T 9l Q). In 
addition, using the canonical identification between the vector spaces Tj Q and fcer(rg)^ =< 
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d d 
—r > (respectively, T qi Q and fcer(rg)^ 1 =< — — >) we have that the 1-form D 2 L(v ,vi) 

(respectively, D±L{vo,v\)) may be considered as the 1-form on Q at the point go (respectively, 

i=l a ^0|(« o ,t;i) 

(respectively, ^ ^ dq\{ qi )). 

i=l ^1 |(«Oi"l) 

5.1. Variational calculus on Cjv,[o,T]- Define the semi-discrete action § S( ^L : Cjv,[o,t] — * K 
as follows 



dt 



Jo T,k=o L ivi^ *)> *)> y( fc + 1. + 1, *)) 



dt 



Definition 5.1. An element y € Cjv,[o,T] * s a solution of the semi- discrete variational problem 
determined by L if and only if it is a critical point of the Lagrangian system defined by § s ^i- 

Therefore, a solution y of the semi-discrete variational problem extremizes S sc ]£ among all the 
possible variations of y, where a variation of y is a smooth curve s £ (— e, e) — ► y s G Gn,[o,t] 
with y — y. Denote by 

Sy k (t) = ^(k,t) 



ds 



s=0 



where we use the notation y k (t) = y(k,t). Observe that y s (0,t) — yo(t) and y s (N,t) = j/jv(i) for 
all s. We will also use the following notation 

Ai(i,i+i)(t) = {D l L){y k {t),y k {t),y k +i{t),yk+i{t)) , 

DjL {ktk+l) (t) = (DjL)(y k (t) , y k (t) , y k+1 (t) , y k+1 (tj) 

forte {1,..., 4} and je {1,2}. 

Extremizing the semi-discrete action function among all the possible variations, we find that 



d_ 

ds 



s =o ^L(y s )) = Ts 



s=0 



N-l 



E L {Vs{k,t),y s {k,t),y s {k + l,t),y s (k + l,t)) 



k=0 



dt 



T N-l 



D 3 L( k -i, k )(t)5y k (t) + D 4 L( k _i ik )(t)—5y k (t) 



+D 1 L {ktk+1) (t)5y k (t) + D 2 L {ktk+1) (t)—8y k {t) 

r T [N-l 



dt 



f i d \ 

= J S (^(fe-i.kjW + DiL ( fc,fc+i)(t) - j f {D 4 L {k _ ltk) (t) + D 2 L {ktk+1) (t))J 5y k (t) 

N-l rp 

+ {DiL (k _ hk) (t) + D 2 L {kik+1) (t)) Sy k (t) . 



dt 



k=l 



Therefore, the serai- discrete Euler- Lagrange equations are: 

D 3 L (k _ hk) (t) + I>iL (fcifc+1) (t) - j t (D 4 L (fc _i ifc ) (t) + D 2 L (kik+1) (t)) = 0, 1 < k < N - 1 

(5.1) £>4i(fc-i,fc)(0) + D 2 i(fc,fc+i)( ) = 

£ , 4i(fe-i,fc)(7 1 ) + D 2 L( kk+1) {T) = 

The first equations represent a system of second order differential equations of the form: 
Fiyk-xitf^ki^^k+xi^iyk^x^^ki^^k+xtt);^-!^),^^^),])^!^)) = 
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When the matrix {D^iLu-ij) (t)) is regular then we may locally write these equations as 
(5-2) Vk+i(t) = G(y fc _i(f),y fc (t),y fc+ i(i);y fe _i(f),y fc (t),y fc+ i(t);y fe _i(f), y k (t)) 

Then, for enough small T, the semi-discrete flow T: 

T : Q[°< T 1 x — » x 

(yk-i(-),Vk(-)) 1 — ► {vk(-),Vk+i(-)) 

is well-defined (since Equations (|5.2p appear as a system of explicit differential equations) , where 
(yk-i{-),Uk(-),yk+i(-)) satisfies Equations (fSTTj) and 

Q[°< T 1 = {z : [0,T] Q/z e C 2 ([0,T])}. 

5.1.1. Symplecticity. Define the Poincare-Cartan 1-forms 9^, 9j £ A 1 (<5'°' T ' x Q[° ,T 1) as follows 
9£(X(.),y(0) = -Jo(DiL (OA) {t)X(t))dt 

= - / T (DjI^dWIW + D 2 L (0 ,i)(t)X(t)) A 

- / T (£>i£ ( o,i) (*) - I (£> 2 £ (0 ,i)(i))) dt - £> 2 L (0il) (t)X(t) 

e+(x(-),y(-)) = (D 2 L {0A) (t)Y(t)) dt 



So [D 3 L {QA) (t)Y(t) + D A L {0tl) (t)Y(t)j dt 

T 




Jo (D 3 L {QA) (t) - | (D 4 L {0 ,i)(t))) r(# + D 4 £ (0il) (t)y(t) 



Then, we have 

e+(x(.),r(-))-ez(x(.),y(-)) = d 



T 

£(ifo(*),w>(*),ift (*),»!(*))# 



Therefore, there exists a well-defined 2-form fix, 

fi£ = d®l = d&t 
As a consequence of Equations (|5.1|) we deduce that T*9£ = 6 \ and 

5.1.2. Legendre transformations. Define the semi-discrete Legendre transformations as 

Legl : Ql°' T l x Q[°- t 1 ► t*q[o,t] 

(2/o(-),2/i(0) — Le^(yo(-),yi(-)): ^ (.)Q [ °' T1 — » 

A(-) 9£(X(-),0) 

and 

Legl '■ Q l °' T] x Q l "' T] * T*Q [0,TI 

y(-) —> e+(o,y(.)) 

Denote by 9q and 17q the Liouville 1-form and the canonical symplectic 2-form, respectively, on 
T*Q[°< T ] defined by 

6q(«(.))(1(-)) = 5(-)((r Q[ o. rl )*W-)))=/o T «((^)*(^W))* 

= d&Q 

where tq : T*Q — > Q and Tq[o,t] : T*Qt°' T l — > Q[°> T ] are the canonical projections. Then, it is 
easy to prove that 

{Leg- L )*Q Q = @l , (Le<^)*9 Q = 6+ (Le<^)*Q Q = fi L . 
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5.1.3. Momentum mapping. Suppose that the Lagrangian is invariant by a Lie group of symme- 
tries, that is, if cf) : G x Q — > Q is the action of a Lie group then 

L{T q (f) g (v q ),T ql (j> g {v q t)) = L(v q ,v q /), Vv q , v q > G TQ, q,q'eQ 

Infinitesimally, this condition means that 

ZT(QxQ)(v q ,Vg<)(L) = (£ T qK) +Ctq(v)) ( L ) = °> V £ e 9 

where £tq is the infinitesimal generator of the lifted action <fi T : G x TQ — ► TQ. In particular, 

KTQ(i/fc-i(*),yfc-i(t)) + tT Q (yk(t), y k (t))) (L) = 

or, in other words, 

0i£(fc-i,k)(*)6j(lfo-i(*)) + D 2 L {k _ hk) (t)- (e Q (» fc _i(t))) 

+£>3i(fc-i,fc)(*)€g(tffc(*)) + DtL {k _ lM (t)j t (£ Q (y fc (i))) = 0, 

where £q is the infinitesimal generator of the action </> corresponding to £ (note that £tq is the 
complete lift of £q). Therefore 



£>3i(fc-i,fe) (*)£q (l/fc (*)) + (*)&} (yk-i (*)) 

-- (D 4 L (fc _ lifc) (t)) &}(»*(«)) - -| (I?2£(fc_i,fc)(*)) fg(l/fc-i(*)) 

+^ [D4£(fc_i lfc )(t)eo(i/fc(i)) + £> 2 i (fc _i, fc )(i)eg(|/fc-i(t))] = 
Subtracting first expression in (|5.ip applied to £t(Qx<2) and the above equation we deduce that 

DiL (k<k+1) (t)Z Q (y k (t)) - j t (D 2 L [Kk+1) {t)) Z Q (y k (t)) 

= D l L {k _ hk) {t)ti Q {y k _ 1 (t)) - j t (D 2 L {k _ hk) {t)) f Q (tffc-i(t)) 

+ ^ [£>4i(fc-i,*)(t)^(l/fc(t))+i>a£(fc-i,*)(t)^(tffc-i(t))] • 

Therefore, integrating and using the two last equations in (|5.ip we obtain the following preser- 
vation law 

rT r j 



£>l£(fc-l,fc)(*) - jT (-D2i(fc-l,fc)(*)) 



foOte-iW) dt + D 2 L {k _ ltk) {t)i Q {y k ^ l {t)) 



d 



DiL( k k+l )(t) - — (L) 2 i(fe ! fc+i)(i)) 



i Q {y k {t))dt + D 2 L (k , k+l) (t)t Q {y k {t)) 



Note that this equation may be written as 

Z(2/fc-i(-):yfe(-))(CQ °yk-i,£,Q oyk) = <dI(yk(-),yk+i(-))((,Q °yk,£,Q °Vk+\)- 



5.1.4. The nonlinear wave equation. First point of view. The semi-discrete Euler-Lagrange equa- 
tions for the Lagrangian: 

2 



L sd (u ,ui, (tt )t, (ui)t) 



1 / («o)t + («i)t 



(itfc-i)tt + 2(it fc )tt + (Mfc+i)tt _ f_ 
4 ft, 



, / ~ Ufc \ i I u k — U k -i 

— J"" v ft 

0, 1 < Jfe <N-1 
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with boundary conditions 

(5.3) 
(5.4) 



Uo(t) and u±(t) fixed 

(«fc-i)t(0) + 2(«fc)*(0) + (« fc +i)t(0) = 

(« fc -i)t(T) + 2(u fc ) t (T) + (ufc+i) t (T) = 



5.2. Variational calculus on Sjvjo.Tl- Instead of the set of functions Gn,[o,t] we consider the 
set Sjv,[o,t] where now the values at time and T are fixed. As in the previous case, we have the 
semi-discrete action S sc ji : Cjv,[o,T] — > ^: 



7V-1 



L(v(k, t), y(k, t), y(k + 1, t),y(k + 1, t)) 



fc=0 



dt . 



Definition 5.2. A function y G Cjv,[o,T] * s a solution of the semi-discrete variational problem 
determined by L if and only if it is a critical point of the Lagrangian system defined by S SC [L. 



Therefore 
d 



ds 



(Wys)) = 



d_ 

ds 



N-l 



^2 L iys{k,t),y s {k,t),y s (k + l,t),y s (k + l,t)) 



k=0 



dt 



T N-l 

is 



D 3 L {k _ l!k) (t)5yk(t) + D 4 L (k _ hk) (t) — 6y k (t) 



+DiL( ktk+1) (t)8y k (t) + D 2 L( ktk+1) (t)—5y k (t) 



+ 



+ 



DiL ( o,i)(t)tflto(t) + D 2 L (0A) (t)-Sy (t) 



dt 
dt 

dt 



D3L( N -i, N )(t)5y N (t) + D 4 L {N _ 1 ^ N) (t)—5y N (t) 



f 

Jo 

-T VN-1 / d \ 

(fe-i,fe)(*) + D 2 L^ k k+1 - ) (t)) ) 8y k (t) 



dt 



t r 



DiL {0A) (t) 



dt 



DoL 



2-^(0,1) 



(*)) 



dya (t) dt 



+ 



£>3-£(iV-l,iV)(i) - (At-k(jV-l,iV)(i)) 



Therefore, the semi-discrete Euler- Lagrange equations are: 

^3-k(k-l,fc)(*) + ^l-t-(fc,fe+l)(*) - ^ (^^(fe-l.fe) (i) + D2i(fc,fc+1)(<)) 



(5.5) 



D 1 L {0A) (t)--(D 2 L {0 , 1) (t)) 



^3-^(iV-l,JV)(*) 



fit 



(N-1,N) 



(*)) = 



0, 1 < fc < JV-1 







These are precisely the Euler-Lagrange equations 

&L d { dL^ 



0, < k < N 



dy k dt \ydy k ^ 
for the Lagrangian Z : (TQ) Ar+1 ~ TQ N+1 — ► R defined by 



N-l 



L(y , ...,y N , (yo)t, (yjv)t) = X! L ^ fc ' (Vk)t,yk+i, (yk+i)t 



k=0 
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where (TQ) N+1 stands for the Cartesian Product of N + 1 copies of TQ. 
We define 

• The Poincare Cartan 1-form 6 Z G A 1 (TQ N+1 ): 

N-l 

&l = D 2 L (0A) dy Q + ^2 [ D i L {k-i,k) + D 2 L {k ^ k+1) ] dy k + D 4 L (N _ ljN) dy N 
fe=i 

• The Poincare- Cartan 2- form f2^ = g?0£. 

• The energy function 

AT-l 

E~ L = £> 2 i(o,i)(yo)i + 51 [-D4i(fe-i,fc) + £>2-k(fe,fc+i)] (yfc)i + D i L^ N _i, N ){y N ) t - L 

k=l 

We say that the system is regular if and only if the 2-form fi^ is a symplectic 2-form. Locally the 

( d 2 L \ 

regularity is equivalent to the non-singularity of the Hessian matrix — — r— —. — — < k, I < N . 

\d{yk)td{yi)t ) 

In such a case, there exists a unique vector field ^ E X(T(Q) N+1 ) such that 

(5.6) it L n L =dE L 

Moreover, the integral curves of £g are the tangent lifts of the solutions of the semi-discrete Eulcr 
Lagrange equations (|5.5[) . 

In many situations, the 2-form f2^ is not symplectic. Then, (|5.6[) has no solution, in general, 
and even if it exists it will not be unique. Let bi : TQ N+1 — > T*Q N+1 be the map defined by 
bi(X) = ix^i- It may happen that is not surjective. We denote by kerO^ the kernel of 6^, 
i.e., kerfo^ = kerfJ^. 

In [H [5] , Gotay and Nester have developed a constraint algorithm for presymplectic systems 
which is an intrinsic version of the classical Dirac-Bergmann algorithm . They consider the set of 
points P 2 of Pi — TQ N+1 where (|5.6p has a solution and suppose that this set P 2 is a submanifold 
of P\. Nevertheless, these solutions on P 2 may not be tangent to P 2 - Then, we have to restrict 
P 2 to a submanifold where the solutions of (|5.6[) are tangent to P 2 - Proceeding further we obtain 
a sequence of submanifolds: 

> P k >P 2 ^Pi= TQ N+1 . 

Alternatively, these constraint submanifolds may be described as follows: 
Pi = {P € Pt-i I dE L (p)(v) = , to G T p PU } , 

where 

T P Pti = T pPi I ^l( x )( u ' v ) = , Vu e TpP^ } . 
We call P 2 the secondary constraint submanifold, P3 the tertiary constraint submani- 
fold, and, in general, Pi is the i-ary constraint submanifold. 

If the algorithm stabilizes, i.e. there exists a positive integer k G N such that P k — P k +i and 
dimPfe 7^ 0, then we have a final constraint submanifold Pf = P k , on which exists a vector 
field X such that 

(5.7) (i x n L = dE~ L ) /Pf . 

If £ is a solution of (|5.7p then every arbitrary solution on Pf is of the form £' = £ + Y, where 
Y G (kerfi^n TP } ). 

Another interesting aspect of this theory is that, in any case, regular or singular, since the 
Lagrangian L is autonomous, then £{E£) = for any solution of Equation (|5.7p . In particular, 
for we have that 

N-l 

•D2i(oa)2)o + X! [ D * L {k-i,k) + £>2£(fc,fc+i)] Vk + D 4 L {N _ 1 ^y N - L =0 



df 



k=l 
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and the energy is a constant of the motion. 

5.2.1. Legendre transformation. Define the semi-discrete Legendre transformation as 

FL : TQ N+1 — ► T*Q N+1 

(y a ,...,y N ,(y Q )t,...,(yN)t) ' — > (yo, ■ ■ ■ , Vn;Po, ■ ■ ■ ,Pn) 

where 

Pa = ^2-^(0,1)) ■■•)Pk = D4,L(k-i,k) + D 2 L( k k+1 ), ...,p N = At£(jv-i,iV) 

If Aqjv+i and loqn+i are the Liouville 1-form and the canonical symplectic 2-form on T*Q N+1 
then 

FL*A Q «+i = (-), . ¥L*uj q n + i = <>, 

Moreover, if the Legendre transformation FL is a global difFeomorphism, we will say that the 
Lagrangian L is hyperregular. In this case it is well defined the function H = o (FL) -1 . 

Therefore, we have a Hamiltonian representation of the equations of motion (|5.5[) : 

ix„^Q«+i = dH 

or, in coordinates, 

9H . dH n 
Vi = Pi = -tt-, < i < N 

opt oy % 

Obviously, and Xh are FL-related, i.e. (FL)»^ = Xh- 

If the lagrangian L is singular, FL is not a difFeomorphism. However, we may assume that L is 
almost regular, i.e., M\ — ¥L(TQ N+1 ) is a submanifold of T*Q N+1 and, FL is a submersion 
onto Mi with connected fibers. The submanifold Mi will be called the primary constraint 
submanifold. 

Since the Lagrangian is almost regular, the energy E^ is constant along the fibers of FL. 
Therefore, E^ projects onto a function Hmi on M±: 

H Ml (¥L(pj) - E L (p) ,Vp E TQ N+1 . 

If we denote by i Mt : Mi — > T*Q N+1 the embedding of Mi into T*Q N+1 , then we obtain a 
presymplectic system (Mi, (jmJ'wqw+i, dH^). If we now apply the constraint algorithm to it, 
we shall obtain the following sequence of constraint submanifolds: 

► M k -» • .W_, • Mi , 

as in the Lagrangian side. Denote by Mf the final constraint submanifold (if it exists) for this 
presymplectic system. In Mf there exists at least a vector field X g X(M/) such that 

(ix [(i Ml )*VQ*r+i] = dH Mx) /Mf 

The Legendre map restricts to each submanifold Pj, i > 1, of TQ N+1 and then we obtain a 
family of surjective submersions (FL), : Pi — > Mj which relates the constraint submanifolds Pf 
and M^ in particular, Pj and M/. 

5.2.2. Momentum preservation. Suppose, as in Subsection l5.1.3l that the Lagrangian L is invari- 
ant by a Lie group of symmetries. Infinitcsimally, this conditions implies that 

£TQ« + i(%>, ■ • ■ ,Vq N )(L) = 0, V^eg 

In such a case, applying the classical Noether theorem, we deduce that the function F : TQ N+1 — ► 
K defined by 

F = 0£(Ctq«+i) 

is a constant of the motion. Explicitly, 

F(yo, Vo,---, Vn, Vn) = D 2 L {0A) £, Q (y ) 



N-l 

E 

k=l 



[D 4 L( fc _i ifc) + L» 2 L (fc ^Q(yk) + D 4 L^ N _ 1N ^Q(y N ) 
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5.2.3. The nonlinear wave equation. Second point of view. Consider the semidiscrete Lagrangian 

2 



L ad (uo,ui, (u ) t , (ui) t ) 



1 / (uo)t + (ui)t 



2 V 2 J ' V h 
for the nonlinear wave equation. Then the Euler-Lagrange equations are: 



(«) 1 ((«o)tt + (ui)tt) - ^ 



(b) 



(c) 



(u k -i)tt + 2(ufc) tt + (ufc+i)tt _ 1 
4 /i 



5/'(^) = o 



/ / Wfc+i - Ufe \ r I Uk — Uk-i 
—h—)- a { h 

0, 1< fc< iV-1 



7 ((wjv)** + ("v-i)tt) + -r 

4 a 



(7 



, / U N - UM-1 



The Lagrangian X is 

L(u , . . .,u N ; (u )t, (ujv)t) 

and the Poincare-Cartan 1-form is 

e, . I 



N 

E 

fc=i L 



1 ( (wfc-i)t + (wfc)t 
2 



/ u k - u k --i 
V h 



f 



N-l 



{(uo)t + («i)t) + 2J (( u fc-i)t + 2 ( u fc)« + («fe+i)t) rf "fe 
fc=i 

+ (K-i)i + (u N ) t ) du N ^ 
In this case the Poincare Cartan 2-form ft^ = — dQ>i is degenerate with 



N 



. d 



N 



kerfi £ = span <( — ,£(-1 



. fc=o 



fc=0 



d{u k )t 



and, therefore, it is necessary to implement the constraint algorithm. We will study in a future 
paper the constraint algorithm for this and other examples. 



6. Conclusions and outlook 

In this paper we have elucidated the geometrical framework for semi-discrete Mechanics, an 
useful tool for numerical simulation of Classical Field theories. For simplicity, we only have 
studied the case of the lagrangian L s d : TQ x TQ ~ T(Q x Q) — > R but also it is possible 
to consider extensions of this theory to the case of Lagrangians L s d ■ 3T G — > R where 7 T G 
is the prolongation of a Lie groupoid G =4 M over the fibration r : AG — > R, AG being the 
Lie algebroid associated to the Lie groupoid G (see [5] for more details). This extension may be 
useful for generating numerical schemes for Field Theories modeled on Lie algebroids (see |14j). 

Moreover, the semi-discrete lagrangian L s d : TQ x TQ — ► R is adapted for field theories 
whose lagrangian L : Z — > R does not depend on the base variables (the base space is assumed 
bi-dimensional) , that is, L — L(y l ,z^). Suppose, for simplicity, that the parameter space is 
X = R x P, with P = R 2 and Y = X x R, and the continuous Lagrangian is of the form 



L — L{t,x ,x ,y, Zt, z Xl , z X2 ) 
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Introduce the following natural semi-discretization of L 

(A><i)(ki,fa)(*>y(ki,fc2)(*)>#(fci+i^^ 

tu i u i u y^iM)^) + J/(fei+i,fc 2 )(i) + V{k 1 ,k 2 +i){t) 

— L(t, K1A1, K 2 n-2, g , 

V(k u k^){t) + y(k 1 +iM)( t ) + y(fci,fc 2 +i)(0 2/(fci+i,fc 2 )W ~ y(k u k 2 ){t) yj^M+i)^) - y(fci,/c 2 )ft) x 

3 fti ' /i 2 J 

based on a triangularization of the space Z x Z. Then we need to extremize a functional of the 
type 

r T 

&sd L = 



( i sd)(ife 1 ,fc 3 )(*,2/(fci,fc 2 ) +i,fc 2 )(*)'f(A! 1 ,fc2+i)(*); 

(fc 1 ,fc 2 )6{0,...,iV-l}x{0,...,JV-l} 



2/(fci,fc 2 ) (*)) y(fei+i,fe 2 ) 2/(fei ,fe 2 +i) (*)) 



Observe that (-£ad)(fc 1> k 2 ) g C 2 (1R x TIR 3 ) and then new tools must be used as, for instance, 
cosymplectic geometry and suitable adaptations of the results about Discrete Field Theories (see, 
for instance, [131 118j ). Also extensions for the case of a non-uniform grid may be considered and 
in a future paper will be discussed. The generalization of this theory to the case of constrained 
multisymplectic field theories, as for instance in the case of incompressibility constraints in fluids 
, may be obtained using and adaptation of constrained discrete variational calculus (see pQ , 
for instance). 
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